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Abstract 

We present a calculation of interference effects in Hjj production via gluon fusion 
and via vector boson fusion, respectively, beyond tree level. We reproduce results re- 
cently discussed in the literature, but go beyond this calculation by including a class 
of diagrams not considered previously. Special care is taken in developing a numeri- 
cally stable and flexible parton level Monte-Carlo program which allows us to study 
cross sections and kinematic distributions within experimentally relevant selection cuts. 
Loop-induced interference contributions are found to exhibit kinematical distributions 
different in shape from vector boson fusion. Due to the small interference cross section 
and cancelation among different quark flavor contributions their impact on the signal 
process is found to be negligible in all regions of phase space, however. 
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1 Introduction 



Higgs production via weak boson fusion (WBF), i.e., the reaction pp^Hjj, mediated by 
t-channel weak boson exchange, constitutes a particularly promising production mechanism 
for the Higgs boson. Due to the distinctive signature of two hard jets accompanying the 
decay products of the Higgs boson, this channel is discussed as possible discovery mode for 
a scalar, CP-even boson as predicted by the Standard Model (SM) [H [2], and as powerful 
tool for a later determination of its couplings [3j. Furthermore, WBF could be employed in 
studying deviations from the SM expectations and help to spot signatures of physics beyond 
the standard model. This can only be achieved, however, if accurate measurements are 
matched by precision calculations of SM signal and background processes and predictions 
for possible new physics scenarios. 

Next-to- leading order (NLO) QCD corrections to the SM pp^Hjj WBF signal are avail- 
able for cross sections [I] and distributions At the same level of accuracy, some of the 
most important background processes such as Vjj [6] and VVjj [7] production in WBF, 
tt [8] and ttj [9] production are known. Beyond the standard model, Monte Carlo studies 
have been performed for WBF Hjj production in the presence of anomalous gauge boson 
couplings [TO] and in the context of supersymmetric models [TT]. Recently, NLO electroweak 
(EW) corrections for cross sections and distributions have been presented [12] . Finite parts 
of the NNLO-QCD corrections have been calculated in [13] and found to be negligible. 

An irreducible background to the Higgs signal in WBF is constituted by Hjj production 
via gluon fusion (GF). Higgs production via GF is mediated by a heavy quark loop. An exact 
calculation of pp—>Hjj via GF at the lowest non- vanishing order has been performed in [TT] . 
NLO-QCD calculations have made use of the large top mass limit, where the coupling of 
the Higgs to gluons is parameterized by an effective vertex [15]. Phenomenological studies 
have revealed the complementary features of the WBF and the GF Hjj production processes, 
suggesting search strategies for suppressing the GF channel as background to the clean WBF 
signature. More recently, GF has also been considered as a signal process [TB], because the 
dijet angular correlation is sensitive to the CP parity of the Higgs boson. 

Although GF and WBF are usually considered as separate reactions, their interference 
in the qq^qqH subprocesses is possible. In Refs. [TTJ [18] interference effects at tree-level 
due to identical flavor effects have explicitly been shown to be tiny and entirely negligible 
for cross sections and distributions. The authors of [IB] speculated that loop-induced inter- 
ference effects should be large. An explicit calculation revealed, however, that loop-induced 
interference effects are also small [19J. 

We present a similar calculation for the same process, studying interference effects be- 
tween GF and WBF which emerge beyond tree level. Being performed with entirely different 
methods, our work confirms the main findings of Ref. [IS] and extends it in two respects: 
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Figure 1: Feynman diagram for the tree level process qq'^qq'H via WBF. 

First, we include a class of real emission contributions which has been neglected in [19]. Sec- 
ond, we develop a fully flexible parton level Monte Carlo program which allows us to study 
cross sections as well as arbitrary distributions within experimentally feasible selection cuts. 
Being implemented in the modular vbfnlo environment [20 J, the impact of the interference 
contributions on the WBF signal is studied in detail. 

We give a thorough outline of our calculation in Sec. [21 Numerical results are discussed 
in Sec. El We conclude with a brief summary in Sec. HI 

2 Elements of the calculation 
2.1 General framework 

Hjj production in WBF mainly proceeds via quark scattering, qq^qqH . In Ref. [T7| con- 
tributions from identical-flavor annihilation processes such as qq-^Z* —>ZH with subsequent 
decay Z^qq or similar WH production channels have been shown to be entirely negligible 
in the phase-space regions where WBF can be observed experimentally. In the same work, 
identical quark interference effects from qq^qqH and crossing-related channels were demon- 
strated to affect cross sections and kinematic distributions at an insignificant level. This 
finding was confirmed by Ref. [TS]. In the following we will therefore restrict our discussion 
to quark scattering via exchange of a weak boson in the t-channel, i.e. the reaction qq'-^qq'H, 
where q and q' stand for quarks of different flavor, see Fig. HJ We will refer to the respective 
tree level scattering amplitude by A^^ BF . Color factors are not included in -M^/bf an< ^ 
will be denoted separately. Adaptation for the crossed processes qq'^qq'H, qq'^qq'H, and 
qq'—>-qq'H is straightforward. 

Higgs production in quark scattering reactions, mediated by a gluon in the i-channel 
which couples to the Higgs boson via a top-quark loop is depicted in Fig. 12(a). For a Higgs 
mass well below the top-pair production threshold, the coupling of the gluon to a scalar, 
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Figure 2: Feynman diagrams for the tree level process qq'-^qq'H via GF, mediated (a) by a 
top-quark loop and (b) by the effective Hgg vertex of Eq. ( 12. II) . 



CP-even Higgs boson can be parameterized by an effective Lagrangian of the form 

a 
12-kv 



£eff = ^r.HG%G^ , (2.1) 



where G denotes the gluonic field tensor and v = 246 GeV the vacuum expectation value 
of the Higgs boson. The respective Feynman diagram is depicted in Fig. 12(b). Throughout 
this work we will employ this effective coupling for the Hgg vertex. In the following, we 
denote the lowest order scattering amplitude for qq'^qq'H via GF by .Mgj" Anal ° 

ETOUS tO 

the WBF case, color factors are not included in the amplitude M-qp- 

At tree level, the GF and WBF production channels for qq'—^qq'H do not interfere due to 
the color structure of the two processes. An interference between GF and the neutral- current 
contributions to WBF becomes possible, however, if an additional gluon emission is consid- 
ered. Flavor-changing H^H^-fusion diagrams cannot interfere with the flavor-conserving gluon 
exchange diagrams. For the neutral-current mode, two types of loop contributions emerge: 



(1) One- loop diagrams, where a gluon is exchanged between the upper and the lower 
fermion line in the WBF diagram of Fig. [1] (for V = Z). The respective loop amplitude 
-^wbf° P ' ) yields non-vanishing contributions at order O{a 2 o i s ) when interfering with 
the tree-level GF production amplitude M^q F - Here, we count the HZZ coupling as 
av and the Hgg coupling as a s /v; see Fig. |3]for a representative Feynman graph. 

(2) GF diagrams with an extra gluon exchanged between the upper and the lower quark 
line, A1q F 1oop \ also contribute at order 0(a 2 a z s ) when interfering with the tree-level 
ZZ-fusion amplitude A^bf as depicted in Fig. HI 



All relevant loop diagrams involve pentagon diagrams. Box, triangle, and bubble diagrams 
do not contribute due to color conservation. 

At the same order in the perturbative expansion, real emission diagrams have to be 
considered. Non-vanishing neutral-current contributions to the qq'^qq'gH process arise 
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Figure 3: Representative loop contribution to the interference cross section for qq'—^qq'H via 
WBF and GF, respectively, at order 0(a 2 a^), where the 1-loop WBF amplitude interferes 
with the tree-level GF amplitude. 

q » — -j s> =» c » q 



q' > — 3 ^ > » » q' 

Figure 4: Representative loop contribution to the interference cross section for qq'—^qq'H 
via WBF and GF, respectively, at order 0(a 2 a^), where the 1-loop GF amplitude interferes 
with the tree-level WBF amplitude. 

from the interference of scattering diagrams like those depicted in Fig. [51 The upper diagram 
shows an interference between gluon emission from the g-line in the WBF amplitude and from 
the g'-line in the GF amplitude. In the lower diagram the inverse configuration is illustrated. 
Interference graphs where both gluons are emitted from the same quark line cancel out when 
colors are summed over. The same applies to those graphs of the GF amplitude where a 
gluon is attached to the internal gluon line or to the Hgg vertex. We denote the real emission 
amplitudes for qq'^qq'gH that do not cancel by .M^bf an d respectively. Further 

contributions to qq'^qq'gH scattering via GF, referred to as M^p , arise from a topology 
absent in qq'-^qq'H, where the Higgs boson is radiated off the final-state gluon rather than 
the t-channel exchange boson, see Fig. [HI 

Since only diagrams with gluons being emitted from different quark lines contribute to the 
real emission, no collinearly divergent configurations emerge. Singularities arise, however, 
when the final-state gluon in Ai^BF or -A^p is soft. Such divergences in the real emission 
diagrams are eventually canceled by respective singularities in the virtual contributions. To 
isolate them in intermediate steps of the calculation, a proper regularization scheme has to be 
utilized. We therefore perform our calculation in d = 4 — 2e dimensions, implementing both 
the dimensional regularization and the dimensional reduction prescriptions [21j. Checking 
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Figure 5: Representative cut amplitudes for the qq'^qq'gH process via the interference of 
WBF and GF amplitudes, at order 0(a 2 a 3 s ). 




Figure 6: Representative diagram contributing to the interference of qq'—^qq'gH via WBF 
and GF, respectively, at order 0(a 2 a 3 s ). 

that both schemes yield the same results provides a test of our calculation. 
2.2 Virtual contributions 

For the discussion of the loop contributions to Hjj production via WBF and GF, respectively, 
we resort to the quark-quark scattering process 

q(Pa) + q'(p b )^q(pi) + q\P2) + H(p H ) . (2.2) 
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The one-loop amplitudes we are considering, 
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-^WBF OP) = E^ WBF (Pa,Pi;p 6 ,P2)^ fc , (2.3) 



k=l 
13 



M (i-ioo P ) = Y,F^(p a , Pl]Pb ,p 2 )M k , (2.4) 



k=\ 



can be expressed as linear combinations of process-dependent pre-factors, .F^bf an( j ^gf^ 
and fermion spinor chains Xik, so-called "standard-matrix elements" (SME). Following 
Ref. [22], we introduce 



F {<U/? 7 } = U (Pli X l) {la,7a7f3l~,}u(p a ,\ a ) , (2.5) 
r {'i/37} = «(P2, A 2 ) {7a,7a7/37 7 }^(P^ ^ ( 2 - 6 ) 

where u(pi, Aj) denotes the quark spinor for fermion % with momentum p^ and helicity Aj = 
±1/2. For contractions with an arbitrary momentum p we use the shorthand notation 
T p = TuP 11 . For the reaction (12.21) . 13 SME emerge, 

M {1 ,2 } = rw^-^}, m {3A} = r% apb r**>{«*piP.} i 

M {5M = rwr^-wi}, tvi {7i8} = i^iW'W^ 

-^{9,10} = r^r^/^M .m {11 , 12} = r^ P2 r^,WMM 

-^13 = r^ 7 r^',^7. (2.7) 

The SME are computed in two independent ways by means of the helicity amplitude for- 
malism of Ref. [23] and the Weyl-van der Waerden formalism of Ref. [21], respectively. 

The coefficients F^ BF and Ff? F contain coupling factors and remnants of scalar and 
tensor loop integrals, up to rank two and up to five propagator denominators. For the 
computation of the tensor coefficients we have employed two different methods and developed 
two completely independent computer codes. These codes agree with a relative accuracy 
better than 10 -8 for non-exceptional phase-space points away from the zeroes of the Gram 
determinant. The basic features of the two implementations are described in the following. 

Passarino-Veltman type tensor reduction: 

In one of our implementations, we have used the conventional Passarino-Veltman (PV) 
reduction formalism [251 126] for the computation of tensor integrals up to boxes, generalizing 
the method in a straightforward manner to pentagons in the framework of dimensional 
regularization. All tensor coefficients are expressed in terms of scalar master integrals in (4 — 
2e) dimensions with a regularization scale /i. Singularities are manifested as single or double 
poles in e. For GF, only integrals with vanishing internal masses emerge. Expressions for the 
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respective four-point integrals with up to two external off-shell legs are taken from Ref. 
The infrared divergent two- and three-point integrals can be extracted from Refs. [26j l28]. For 
WBF, scalar integrals with up to two massive propagators are needed. We have calculated 
the divergent box integrals by extracting corresponding expressions in mass regularization 
from Refs. [291 30] . Employing the method described in Ref. [28] these expressions are 
then transformed to dimensional regularization (see App. |A} 0. The remaining finite scalar 
integrals were calculated with LoopTools [31], and have been compared numerically to the 
expressions given in Ref. [19j . 

The double and single pole terms of the coefficients F^ BF and are calculated an- 
alytically with the help of mathematica. To this end, we perform the tensor reduction in 
two steps: first for the singular pieces in d = 4 — 2e dimensions, and then in four dimensions 
for the non-singular terms. For tensor integrals up to rank two this separation is trivial, 
since poles do not mix with the finite terms when the reduction formalism is applied. Due 
to the absence of collinear configurations in the interference process we are focusing on, all 
double poles present in divergent scalar integrals cancel in the full loop amplitudes. This 
cancelation occurs for each coefficient Ff? F , F^ BF (k = 1, . . . , 13) separately. The remaining 
single poles for the GF and WBF contributions are proportional to the respective tree-level 
amplitudes such that the singular parts of the loop-induced interference contribution take 
the form 



E 2Re 



..(l-loop) 4.(0)* .v(l-loop) M (0)± 
' WBF ^'GF + JVI GF J * WBF 



smg 



|r(l + e) ^i^2Re 



- /W WBF- /W GF + - /W GF- /W WBF 



fcl J!* -In ^ -In Win' * 12 



AlTfl 2 J \4lTfl 2 ) \4lTfl 2 J \4lTfi 



(2.8) 



where ^2 denotes averaging over initial-state spin degrees of freedom and summation over 
final-state ones. Here, we have introduced the notation s^- = Ipi ■ pj. We will see below that 
the divergent pieces are canceled exactly by respective poles in the real emission contribu- 
tions. In our Monte-Carlo program we thus retain only the finite parts of the scalar integrals. 
The tensor reduction for the non-singular terms can then be performed numerically. 

Denner-Dittmaier type tensor reduction: 

In an independent implementation, we have performed the reduction from five-point to four- 
point integrals by means of the Denner-Dittmaier (DD) reduction formalism [32], while 
still reducing the four-point integrals with the conventional PV tensor reduction discussed 
above. Like the PV method the DD reduction is formulated in d dimensions, but has the 
advantage of avoiding subtle cancelations, which can spoil the numerical integration. Such 
cancelations originate from terms with a very small Gram determinant in the denominator. 



1 We identified a misprint in one of the respective expressions of Ref. [19] 



8 



If two or more external momenta are linearly dependent, the Gram determinant vanishes 
and we encounter numerical instabilities in the respective phase-space regions. In the PV 
approach, one power of the Gram determinant appears, in general, in each step of the 
iterative reduction to lower-rank tensors, and also when an A-point function is reduced to 
(N — l)-point integrals. In contrast to the PV reduction, the DD approach avoids Gram 
determinants in the denominator in the reduction from five-point to four-point integrals. In 
addition, it reduces the rank of the emerging four-point tensor integrals by one. A similar 
reduction formalism has also been given in [55] . 

As expected, the DD method turns out to be numerically far more stable than the PV 
reduction. Throughout our numerical studies, we will therefore resort to the DD formalism. 
The PV reduction is used only to test our results. 



2.3 Real emission contributions 



For the real emission contribution the partonic subprocess 

q(Pa) + q'(Pb)^q(pi) + q'(p2) + g(p g ) + H{p H ) (2.9) 
has to be considered. The amplitudes A^^bf aric ^ ^gf = A^Q F al ' t ' ) + A^Q F al ' f ^ are computed 



numerically by means of the helicity amplitude formalism of Ref . [23] . Results for qq'—^-qq'gH, 
qq'^qq'gH, and qq'^qq'gH are obtained analogously. Due to color conservation, crossing- 
related subprocesses with a gluon in the initial state such as gq'-^qqq'H do not contribute 
to the interference cross section. 

To test our implementation we have compared our results for A^gp and A4q^ at the 
amplitude level and for the interference contribution ^22Re M^^Mq^* + M.q^ A^^bf* 
to madgraph [3lj and found complete agreement within the numerical accuracy of our pro- 
gram. 



2.4 Subtraction procedure 



The real emission contributions contain soft divergences which eventually cancel the corre- 
sponding poles in the virtual contributions, cf. Eq. (12. 81) . A convenient method for isolating 
the singularities is the so-called phase-space slicing procedure. It relies on splitting the 
qq'gH phase space into soft and hard regions by a suitable cutoff parameter and performing 
the integration of the real emission contributions in the two regimes separately. To check 
our results we implement two conceptually different slicing methods: The two-cutoff slicing 
method of Ref. [35] which has been developed in the context of mass regularization, and 
the phase-space slicing method of Ref. [36] which utilizes a Lorentz-invariant cutoff and 
dimensional regularization. 
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2.4.1 Lorentz-invariant phase-space slicing 



The slicing method of Ref. [36] divides the phase space of the final-state particles into a 
hard region where all partons can be resolved and an infrared region for soft and collinear 
configurations. In general, special care is necessary to separate soft and collinear regions in 
order to avoid double counting of singular configurations. Since the interference contributions 
of our interest are free of collinear singularities, the formalism can be greatly simplified, 
however. In the case of qq'-^qq'gH the gluon is considered as infrared when 

s ig = 2pi ■ p g < s min , with % = a, b, 1, 2 (2.10) 



for an arbitrarily small cutoff parameter s mrn , where we closely follow the notation of Ref. 
While for collinear configurations only one Sj g is small, the soft region is defined by requiring 
at least two invariants to be smaller than s m i n . The partonic real emission cross section can 
then be decomposed into a soft and a hard part, 

cereal "soft i "hard /n i i \ 

a = <y + o l^-llj 

The integration over the gluonic degrees of freedom is performed analytically in <3" soft , but 
purely numerically in a hard . 

In order to calculate <r soft we use the factorization properties of the real emission ampli- 
tude in the soft limit. As the energy of the emitted gluon becomes small, the qq'-^qq'gH 
interference amplitudes can be approximated by the tree-level interference amplitudes mul- 
tiplied by a sum of eikonal terms, 



E 2Re 



A/f ( real ) A/f( real )*J_ A/f ( real ) A/f ( rea1 )* 
• /W WBF- /W GF + M GF yva WBF 



soft 



(4vr« s ) ^ 2Re 



2s ab 2s a2 2s bl 2s 12 



SagSbg S ag S 2 g S bg S\ g Slg&2g 



(2.12) 



The color structure of the soft contribution will be considered below. In the soft region, the 
four-particle qq'gH phase space factorizes into a three-particle qq'H phase space and the soft 
gluon phase space for the respective configuration, 

[d(PS 4 )f oft = d(PS 3 ) d(PS g ) soit ^j,g) , (2.13) 

where d(PSs) contains the flux factor, l/(2s), with s denoting the partonic center-of-mass 
energy squared. For two outgoing partons i and j, d(PS g ) soit (i, j, g) is given by [36] 

soft ( 47r ) e sir 1 

d(PS g ^ (i, jig) 167T^ r(l ) ^^^^^^ dSi g dSj g 9(s m i n ^ig) ^(^min ^jg) ■ 

(2.14) 
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The integration over the soft gluon phase space can be performed for each term in the soft 
interference amplitude Eq. (I2.12p explicitely, using 



filP j d(PS 9 r {t (z,j,g) 



2s 



SigSjg 



lir 2 T(l -e) 



2\ e 



;2.i5) 



The generalization of this expression to cases where one of the partons i, j is incoming rather 
than outgoing is straightforward. The soft part of the real emission cross section then takes 
the form 



a 



soft 



CaC f 



[d{PS 4 )Y oit ^2Re iM^Mff* + M^M 



f( real ) A,f( rea1 )* 
l WBF 



- soft 



^V(l + e)^ I d(PS 3 ) ^2Re 




- /W WBF- /W GF + • /vl GF yvi WBF 



1 

+ 2 



In 2 



Sab 



In 



In 2 



Sab 



hi 



Sa2 



l n ,_i^Uln / 512 



Sa2 



In 2 



Sbl 



+ ln 2 



812 



(2.16) 



where we have included the color factor CaCf/2 = 2. When <r soft is combined with the 
virtual contributions, 

a qq ' H = (T virt + <x soft , (2.17) 

all 1/e poles cancel [cf. Eq. (12.81) ] . The remaining terms are finite and can be integrated 
over the three-particle phase space of the qq'H system and convoluted with the parton 
distributions of the incoming fermions numerically. The resulting three-particle contribution, 
<r qq ' H , depends on the unphysical cutoff parameter s m j n . This dependence cancels, however, 
once <7 qq ' H is combined with the hard part of the real emission cross section, 



a 



full 



^hard _|_ ^-qq'H 



(2.18) 



Checking that the full NLO interference cross section is independent of the cutoff parameter 
therefore provides another important test of our calculation. 



2.4.2 Phase-space slicing with energy-cutoff 

The phase-space slicing method of Ref. [35] in principle requires two cutoff parameters for 
separating finite from collinear and soft divergent regions. For our application, however, no 
collinear singular configurations emerge. Thus, applying a single cutoff on the energy of the 
potentially soft gluon is sufficient. 
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In analogy to the Lorentz-invariant slicing method described previously, the real emission 
contribution can be evaluated numerically in the "hard" region of phase space above the 
cutoff, where it is completely finite. Below the energy cutoff, however, the phase-space 
integration over the gluonic degrees of freedom is performed analytically. Following Refs. 
39J, this "soft" contribution to the partonic cross section can be written as 



a 



soft, E— cut 



&s CaCf 

2~ix 2 



d{PS 3 ) J2 2Re 



A/f (0) M yu >* + A// lu; M 

- /W WBF- /W GF + JVI GF JVI WBF 
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Pgl 



e 3 <ae 2ixE 3 
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2s 
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2s 



«.2 



2s 



61 



2s 



12 



^ + 

SagSbg ^ag^2g SbgS\g ^lg^2g 



(2.19) 



where E 3 is the gluon energy, AE the energy cutoff in the rest frame of the two incoming 
partons and A a mass used as regulator. As above, d(PS 3 ) denotes the qq'H phase space. 
Evaluating the integrals over the gluonic degrees of freedom in Eq. (12.191) and rewriting the 
mass-regulated result in terms of dimensional regularization, the corresponding soft cross 
section is of the form 



a 



soft, E— cut 



a s CaCf 
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d(PS 3 ) ]T2Re 
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Sbl 
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gig 



(2.20) 



where the Ei denote the quark energies in the partonic rest frame. In complete analogy to 
the Lorentz-invariant phase-space slicing, the soft contribution to the partonic cross section 
is combined with the virtual cross section. The resulting sum is then free of soft poles and 
can be evaluated numerically. Upon adding the hard part of the real emission contribution, 
the dependence on the cutoff parameter AE cancels. 



2.4.3 Checks 

We have checked that the total pp^-Hjj interference cross section at the LHC within typical 
WBF cuts (for details, see our standard definition of cuts in Sec. [3]) is independent of the 
cutoff parameter for both phase-space slicing schemes. 

For the Lorentz-invariant slicing method, we have varied s m i n in the range 1 GeV 2 < 
Smin < 10 3 GeV 2 . For smaller cutoff values, large logarithms arise and numerical instabilities 



12 




1.5 



-2 1.0 h 

CD 



0.5 z 



0.0 



10 100 1000 

s min [GeV 2 ] 



0.0 



i i i ■ ■ 1 1 1 1 — i i i ■ ■ ■ ■ i 1 — i i 1 1 1 ■ ■ i 1 — i i i ■ ■ ■ 1 1 r 

(b) 



full 



r qq'H 



.hard 



I i I i I 



u 



10 -5 10 -4 10~ 3 10 -2 10 _1 



AE/Vs 



Figure 7: Dependence of the interference cross section for pp^Hjj production at the LHC 
within standard selection cuts on the cutoff of the Lorentz-invariant phase-space slicing 
method (a) and of the energy-cutoff slicing method (b). Shown are cr hard (blue), <r qq ' H 
(green), and their sum, a (red). 



are to be expected. If, on the other hand, a very large value is chosen for s m - m , the soft 
approximation used for determining cx soft is not applicable anymore. Fig. [7J (a) demonstrates 
that the two contributions cr qq ' H and cr hard individually depend on s m - m , while the sum <r fu11 
is constant in the considered range of the cutoff parameter. 

A very similar pattern arises for the energy-cutoff slicing method, depicted in Fig. [7J (b). 
We have normalized the energy cutoff AE by Vs for this study. 



3 Numerical results 



The cross-section contributions discussed above have been implemented in a fully flexible 
parton level Monte-Carlo program, structured analogous to the vbf nlo code [20] which has 
been developed for the study of WBF-type production processes at the LHC. 

The loop-induced interference contributions for Hjj production via GF and WBF we 
consider are a gauge-invariant sub-class of the full NLO-QCD corrections to the scattering 
process pp^Hjj. For the parton distribution functions of the proton we therefore use the 
CTEQ6M set at NLO [10] with a s (M z ) = 0.118. We set quark masses to zero throughout and 
neglect contributions from external top or bottom quarks. As electroweak input parameters 
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we have chosen mz = 91.188 GeV, mw = 80.419 GeV, and the measured value of Gf = 
1.166 x 10~ 5 / GeV 2 . Thereof, we compute sin 2 and a using LO electroweak relations. 
For reconstructing jets from final-state partons, we use the algorithm [41J with resolution 
parameter R kT = 0.8. 

Since we want to study the impact of the interference contributions on the Higgs signal 
in WBF, we apply cuts that are typical for WBF studies at the LHC. We require at least 
two hard jets with 

p Tj >20GeV, |%|<4.5, (3.1) 

where pxj is the transverse component and yj the rapidity of the (massive) jet momentum 
which is reconstructed as the four- vector sum of massless partons of pseudorapidity \rj\ < 5. 
We refer to the two reconstructed jets of highest transverse momentum as "tagging jets" . 
The Higgs boson decay products, which we generically call "leptons" in the following, are 
required to be located between the two tagging jets and they should be well observable. 
To simulate a generic Higgs decay without specifying a particular channel we generate an 
isotropic Higgs boson decay into two massless particles (which represent 77 or bb final states) 
and require 

Pte > 20 GeV , \r] e \ < 2.5 , AR jt > 0.6 , (3.2) 

where ARj£ denotes the jet-lepton separation in the rapidity-azimuthal angle plane. In 
addition, the leptons need to fall between the rapidities of the two tagging jets 

yj,min ^ *^ yj,max • (^"^) 

Furthermore, we impose large rapidity separation of the two tagging jets, 

^Vjj = \Vji -V1/2I > 4, (3.4) 
and demand that the two tagging jets be located in opposite detector hemispheres, 

Vji x Vj2 < , (3.5) 

with an invariant mass 

Mjj > 600 GeV . (3.6) 



To ensure the reliability of our calculation, we have compared our results to those of 
Ref. [19J and found agreement with their main predictions. Diagrams where the Higgs boson 
is radiated off the final-state gluon rather than the t-channel exchange boson as in Fig. [6] 
have not been considered in |19j . This approximation seems reasonable, as we found that 
contributions from these graphs amount to only about 0.3% of the total interference cross 
section. For individual subprocesses they can be larger, however. For the dd^ddH channel, 
for instance, they yield approximately 5% of the subprocess-cross section. 

In Fig. [8] we show the total cross section cr™ t ts for the interference contribution within the 
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Figure 8: Dependence of the total interference cross section cr™ t ts for Hjj production at the 
LHC on the factorization and renormalization scales for the two different scenarios described 
in the text. The factorization scale /if and the renormalization scale /i r are scaled as m# for 
(a) and as the jets' transverse momenta in (b), cf. Eqs. (13. 7p and (13. 8p . respectively. The 
curves show 0f^ t ts as a function of the scale parameter £ for three different cases: £ r = £f = £ 
(solid red), £f = £ and £ r = 1.0 (dot-dashed blue), £ r = £ and £f = 1.0 (dashed green). 



cuts of Eqs. ( I3.1l) - (l3.6p and for a Higgs mass of ran = 120 GeV. The factorization scale, /if, 
and the renormalization scale, /i r , which enters the strong coupling are chosen as follows: In 
panel (a), we set 

/if = ^ra H , a 3 s (// r ) = a 3 s (£ r ra H ) . (3.7) 

In panel (b), we associate the scale for gluon emission from either quark line with the 
transverse momentum of the corresponding jet by setting 

/"f = , «s(A*r) = a s (£rPri) • a s (^p T2 ) • a s {^m H ) . (3.8) 

Due to the absence of collinear singularities, /if enters only via the parton distribution func- 
tions of the incoming fermions, which are mainly probed at rather large values of Feynman x. 
In this regime, the valence and sea quark distributions depend on the factorization scale only 
mildly. Thus, the variation of af^ s with /if is very small. On the other hand, the interference 
cross section exhibits a pronounced dependence on /i r . Since the loop-induced GFxWBF 
interference in qq'—>qq'H production, ofn t ts , represents the first non- vanishing contribution 
in the perturbative expansion, the renormalization scale enters only via the strong coupling 
constant. Thus, the entire /i r dependence of the interference cross section can be traced back 
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initial-state flavor combination 


<t ts [ab] 




NC: + c)(u + c) + (d + s)(d + s) 


51.4 


72.3 


NC: (u + c)(d + s) 


-49.8 


70.8 


CC: (u + c)(d + s) 


- 


405.7 


NC: (u + c)(u + c) + (d + s)(d + s) 


-3.1 


39.3 


NC: (it + c)(d+ s) + (u + c)(d + s) 


2.2 


43.0 


CC: (u + c){u + c) + (d + s)(d + s) 




230.7 


NC: (u + c)(u + c) + (d + s)(d + s) 


4.0 


5.1 


NC: (u + c)(d + s) 


-3.2 


4.3 


CC: {u + c){d + s) 




25.7 


sum 


1.5 


896.9 



Table 1: Contributions of various neutral current (NC) flavor combinations to a^ s (in ab) 
and cr^gp (in fb), and of the charged current (CC) contributions to WBF. Also shown is 
their sum within the cuts of Eqs. fl3.ip - fl3.6|) . 

to the variation of the (/x r ) coupling factor with the renormalization scale. Reminiscent of 
what has been observed for pure WBF production processes (cf., e.g., Ref. [42]) a dynamical 
scale choice as in Eq. (13.81) . see Fig. [8](b), yields predictions with a somewhat reduced scale 
dependence as compared to the fixed scale option of Eq. (13.71) . shown in Fig. [1(a). 

Compared to the total WBF cross section within typical selection cuts, the interference 
contribution we have calculated is almost negligible in magnitude. In Tab. [T] we list a cuts 
for both, interference and pure WBF cross sections for the various flavor combinations of 
the scattering quarks and antiquarks, setting /if = /i r = ™>h- F° r WBF, we consider neutral 
and charged current subprocesses at 0(a 3 ). No ^-exchange diagrams contribute to the 
interference cross section. Tab. [T] reveals the strong cancelations occurring in a^ s among 
the separate channels. While some contributions, in particular for the qq' subprocesses, are 
sizeable, their sum amounts to 1.5 ab only. We will show below that the subtle cancelation 
between the same and opposite isospin qq, qq, and qq scattering contributions leads to 
unexpected shapes of kinematic distributions in flavor-blind experiments. 

Figure M depicts the shapes of the transverse-momentum distributions for "pure" WBF 
Hjj production and for the WBFxGF interference contribution we have calculated. The 
very hard pt distribution encountered for the interference significantly differs from the shape 
of the WBF curve. The small size of the pi distribution at low momentum transfer is mainly 
due to strong cancelations among the different flavor contributions to the full pp^Hjj 
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Figure 9: Panel (a) shows the normalized transverse momentum distributions for the tagging 
jet with the highest for WBF (dashed blue line) and for the WBFxGF interference 
contribution (solid red line). In panel (b) the sum of all positive (dashed green line) and the 
magnitude of all negative contributions (dashed black line) are shown separately. The solid 
red line gives the sum of all contributions, multiplied by a factor of 10. 

interference cross section, as illustrated by Fig. Wh), where the contributions for same isospin 
and opposite isospin qq + qq + qq scattering, <7?° s and cr^ g , are shown separately. The two 
contributions cancel almost precisely to give the total interference contribution, crf^t+cr^t = 
o~i nt . At high pt ; the cancelation effects are less pronounced. In short, the interference 
contribution has a harder transverse-momentum spectrum than expected, because of a very 
efficient cancelation around PT.tag ~ 100 GeV, where the individual distributions peak. 

This cancelation pattern is reflected by the tagging-jet invariant mass Mjj. For studying 
the corresponding shapes of the pure EW and of the mixed QCD-EW production processes, 
we have switched off the invariant mass cut of Eq. (13.61) . The emerging curves are displayed 
in Fig. [10J While the interference cross section is negative at small values of the dijet- 
invariant mass, it is relatively large at high Mjj. Indeed, is is remarkable to find that the 
interference cross section yields an even harder Mjj distribution than the pure WBF cross 
section does. This behavior is somewhat unexpected if considering the rather soft invariant 
mass distribution of the pure GF Hjj production process which has been reported in the 
literature [H] . The full GF pp^Hjj cross section, however, is dominated by gluon- initiated 
partonic channels such as gg—>ggH and qg^qgH. To the interference cross section, on the 
other hand, only quark (and anti-quark) initiated subprocesses contribute, which tend to give 
larger values of Mjj than gluonic contributions. More importantly, the cancelation effects 
reported above in the context of the tagging-jet transverse-momentum distribution affect 
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Figure 10: Panel (a) shows the normalized tagging-jet invariant mass distributions for WBF 
(dashed blue line) and for the WBFxGF interference contribution (solid red line). Panel (b) 
depicts the sum of all positive contributions, cr^° s (dashed green line), the magnitude of all 
negative contributions , — cr^ g (dashed black line), and their sum o"i nt , multiplied by a factor 
of 10 (solid red line). 

the summation over the various flavor contributions to the dijet- invariant mass distribution 
in a similar manner, thereby giving rise to a broad invariant-mass distribution which is very 
small at low values of Mjj. 

The afore-mentioned cancelations have different effects on the rapidity distribution of 
the third, non-tagged jet with respect to the tagged jet located in the positive-rapidity 
hemisphere, 

Vdis = 3/3 - m ax (y x , y 2 ) , (3.9) 

which is shown in Fig. [TT] for the interference contribution and the "pure" WBF cross section. 
The separation of the lowest pr-jet from the tagged jet located in the negative-rapidity 
hemisphere, — 2/3 + min (y x , y 2 ), would be a mirror copy thereof due to our symmetric selection 
cuts. For generating the distribution, we required a minimum transverse momentum of 
PT3 > 10 GeV for the third jet in addition to the selection cuts of Eqs. ( 13.ip - 113.6p . The 
peak of the distribution at small \ydis\ emphasizes that the "soft" jet prefers being close in 
rapidity to the hard jet in the respective detector hemisphere for both, interference and WBF 
contributions. While in WBF the third jet prefers rapidities larger than the associated tag 
jet, j/diff > 0, for the interference contribution y^e peaks at negative values for the various 
flavor contributions and their sum, indicating that the soft jet is typically located in between 
the two tagged jets. This may indicate that the rapidity gap for a color singlet EW-boson 
exchange may in general be filled by the EW-QCD interference contribution. 
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Ydiff Ydiff 

Figure 11: Panel (a) shows the normalized rapidity-separation distribution of the non-tagged 
jet for WBF (dashed blue line) and for the WBFxGF interference contribution (solid red 
line). Panel (b) depicts the sum of all positive contributions, af°^ (dashed green line), 
the magnitude of all negative contributions, — cr^f (dashed black line), and their sum <7i nt , 
multiplied by a factor of 10 (solid red line). 

4 Summary and conclusions 

In this article we have computed the order 0(a 2 al) interference contributions to the 
Hjj production cross section in pp collisions at the LHC via GF and WBF. Since results for 
the total interference cross section and angular distributions have already been discussed in 
the literature [TH] , we have put special emphasis on technical and phenomenological aspects 
of the calculation which have not been discussed elsewhere. In particular, we have given a 
detailed outline of the methods used for the evaluation of loop contributions, the subtraction 
of singularities present in intermediate steps of the calculation, and of the checks we have 
performed to ensure the reliability of our results. In the real emission contributions we have 
included a finite class of diagrams that has not been considered previously. We found the 
numerical value of these contributions small, however. 

Having implemented the interference amplitudes in a flexible Monte-Carlo program based 
on the vbfnlo framework of Ref. [20J, we are able to provide total cross sections and ar- 
bitrary kinematic distributions within experimentally feasible selection cuts. Considering 
the interference cross section as possible "contamination" of the clean WBF Hjj production 
signature, we have studied the associated contributions within typical WBF cuts with widely 
separated hard tag-jets and compared the shape of some characteristic distributions to those 
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of the respective pure WBF curves. We found that, indeed, the interference contributions 
exhibit features rather different from the WBF signal which are caused by strong cancela- 
tions among the separate flavor channels. However, due to the small size of the interference 
cross section which is found to be in the atto-barn range only, the impact of this contribution 
to both, integrated cross sections and differential distributions, is negligible. 
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Figure 12: Momentum and mass assignments for a general scalar box diagram. 



A Infrared divergent scalar box integrals 



In this appendix, we denote the infrared-divergent box integrals with massive propagators 
which emerge in the calculation of the loop corrections to VBF-induced pentagon diagrams. 
We do not list the other scalar loop integrals encountered in our calculation, since they can 
be found elsewhere (see, e.g., [27 | [191 13T] ). 



We obtained the respective soft and collinear singular box diagrams by extracting ap- 
propriate expressions from Refs. [2T?j I5U] in the limit of small quark masses. According to 
Ref. [28J, a relation between different regularization schemes can be established making use 
of the genuine singularity structure of infrared-divergent triangle integrals. With the help of 
this property, we transformed the divergent four-point integrals from mass regularization to 
dimensional regularization. 

In the following, we refer to a genuine scalar four-point function as depicted in Fig. [T21 



A)(<?i, 92, <?3; m 1 , m 2 , m 3 , m 4 ) 



(27T/X) 



A-d 



d d q 



1 



[q 2 — m\ + i5] [(q + gi) 2 — m\ + i5] 
1 



x 



[(<? + qi + qi) 2 — ml + i8] [(q + qi + q 2 + g 3 ) 2 —m\ + i$\ 



J 4 (Si, S2, S3, S4; S12, S23] Tfl 



m 2 ,m 3 , 



ml) 



(A.1) 



where the qi denote incoming momenta of the external legs and the rrii correspond to the 
masses of the internally propagating particles. The kinematic invariants, Si and sy, are 
related to the external momenta via s« = qf and = (qi + qj) 2 . Overlined quantities are 
defined as s = s + i5, etc. 

In this notation, the collinear divergent box integral with two equal internal and two 
different external mass scales, which is sketched in Fig. [T31 takes the form 
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Figure 13: Momentum and mass assignments for the collinear divergent scalar box diagram 
of Eq. ( 1A.2I) . Unlabeled thin lines correspond to massless particles. 
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ln(a&) — ln(a) — ln(6), 

(s — m 2 ){t — m 2 ) + (s 4 — m 2 )m 2 , 

m 2 H {m 2 — s){m 2 — t) + m 2 [m 2 (s + t — s 4 ) + s(s 4 — 2t)], 
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Figure 14: Momentum and mass assignments for the soft divergent scalar box diagram of 
Eq. ( 1A.4I) . Unlabeled thin lines correspond to massless particles. 



k 2 



m 



m 



1 - 4— + 2i5—(m 2 H - 2m 2 )6(Am 2 - m 2 H ) 



m 



H 



m 



H 



L>i(m 2 



m 4 (m 2 — Si) 



D 2 . . [m 2 — s) (m 2 — sA 
— to 



Dim 2 



2m 2 - (1 + p)m 2 . 2m 2 - (1 + p)m 2 „ 

-I— 1Q 

2m 2 2f3m 2 H 



m 2 — s 4 m 2 (t — S4) 

— 2 — 7 + ld ~T~2 — 7\ 
m z — t [m z — t) 



2 ' 



(A.3) 



The soft divergent box integral with one internal and one external mass scale shown in 
Fig. [H] is given by 



/f(si,0,0,0;s,t;0,m 2 ,0,0) 
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s(t — m 2 ) 



r(l + e) f4nfJ?Y | r(l + e) /4vr/i 2x ' 



m 1 - 



In ( — — =■ ) + In 

m z — t 



m 



m 2 — t 



— In 2 - 2 In - In 



m 



-2Li 9 



Si — t 
m? — t 



2Li 2 



111- / \ m 2 — t 
1 '-In 



m 2 — t 



-In 2 



s 1 — 
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(A.4) 



We would like to note that our results agree with those of Ref. [19] . if /i 2 is replaced by 
fx in all terms of Eq. (A. 19) in [19J. 
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